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ABSTRACT 

4U 1822—37 is a Low Mass X-ray Binary (LMXB) system with an Accretion 
Disk Corona. It is one of the very few LMXBs that shows narrow X-ray echpses 
and smaU amphtude pulsations of the neutron star. The X-ray echpse is an excellent 
reference for measurement of orbital evolution of the binary, and we have obtained 16 
new mid-eclipse time measurements of this source during the last 13 years using X- 
ray observations made with the RXTE-PCA, RXTE- ASM, Swift-XRT, XMM-Newton 
and Chandra observatories. These, along with the earlier known mid-eclipse times have 
been used to accurately determine the timescale for a change in the orbital period of 
4U 1822-37. We have derived an orbital period Porb = 0.23210887(15) d, which is 
changing at the rate of Porb = 1-3(3) x 10"^° d d^^ (at To = MJD 45614). The 
timescale for a change in the orbital period is Porb/Porb of 4.9(1.1) x 10^ yr. We also 
report the detection of 0.59290132(11) s (at Tq = MJD 51975) X-ray pulsations from 
the source with a long term average Pspin of -2.481(4) x 10^^^ s s~^, i.e., a spin-up 
time scale (P spin / Pspin) of 7578(13) yr. In view of these results, we have discussed 
various mechanisms that could be responsible for the orbital evolution in this system. 
Assuming the extreme case of conservative mass transfer, we have found that the 
measured Porb requires a large mass transfer rate of (4.2 — 5.2) x 10~® M© yr~^ 
which together with the spin up rate implies a magnetic field strength in the range of 
(1-3) X 10^ G. Using the long term RXTE-ASM hght curve, we have found that the 
X-ray intensity of the source has decreased over the last 13 years by ~ 40% and there 
are long term fluctuations at time scales of about a year. In addition to the long term 
intensity variation, we have also observed significant variation in the intensity during 
the eclipse. Variation was also seen in the pulse profile, which could be due to changes 
in the accretion geometry. We briefly discuss the implications of these results on our 
understanding of the properties of the neutron star and the accretion geometry in this 
source. 

Key words: accretion, accretion disks - binaries: eclipsing, general - stars: individual 
(4U 1822—37) - stars: neutron ~ X-rays: stars 



1 INTRODUCTION 

Low mass X-ray binaries (LMXBs) are stellar systems that 
consist of a compact object, such as a neutron star or a black 
hole, accreting matter from a companion star by Roche lobe 
overflow. Mass transfer in these systems can be conservative 
or non- conservative. In most systems, the mass of the binary 
system and the total angular momentum remains constant, 
which is the case of conservative mass transfer. However, in 
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some systems, a significant mass loss may occur and in such 
systems the mass transfer is called non-conservative. Mass 
loss from the binary system may occur through various pro- 
cesses such as, irradiative evaporation of the secondary star, 
jet emission from the compact star or emission of a wind 
from the accretion disk (Ruderman et al. 1989). Accretion 
may also be driven by the loss of orbital angular momentum 
through gravitational wave radiation, or magnetic braking 
(Hurley et al. 2002; Rappaport, Verbunt & Joss 1983). The 
orbital period of X-ray binaries is expected to change due to 
redistribution of the angular momentum due to interaction 
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between the components of the binary system. Measurement 
of the rate of change of the orbital period (i.e., orbital period 
derivative, Port) of the binary system is therefore necessary 
in order to understand evolution of compact binary systems. 

In some accretion powered X-ray binaries, repeated 
measurements of the orbital ephemeris have led to an ac- 
curate determination of the orbital period evolution. Pulse 
arrival time delay is one of the several techniques used to 
determine the evolution of the binary orbit. Using this tech- 
nique, the orbital evolution timescales have been accurately 
determined in several X-ray binary systems (Her X-1: Deeter 
et al. 1991; Paul et al. 2004; Staubert, Klochkov & Wilms 
2009, 4U 1538-52: Baykal et al. 2006; Mukherjee et al. 
2006, Cen X-3: Kelley et al. 1983; Paul et al. 2007, LMC 
X-4: Levine et al. 2000, SMC X-1: Levine et al. 1993, Woj- 
dowski et al. 1998). Pulse folding and x^ maximization with 
a varying orbital ephemeris have been successfully applied 
to LMC X-4 (Naik & Paul 2004) and SAX J1808.4-3658 
(Jain et al. 2007) to determine the X-ray pulsations and the 
mid eclipse times. However, measurement of Porb for a non- 
pulsing LMXB requires a stable fiducial point in the X-ray 
light curve. Such a measurement is difficult because most of 
the LMXBs do not exhibit sharp orbital features and even 
the pulsating objects often show variable pulse profiles. 

In the absence of pulses from the compact object, eclips- 
ing binary systems provide a good fiducial timing marker for 
precise determination of the orbital evolution. Parmar et 
al. (1986) and Wolff et al. (2009) applied the X-ray eclipse 
timing technique to determine the orbital evolution in the 
LMXB EXO 0748-676. Eclipse timings were also used to 
establish orbital evolution timescales in 4U 1822—37 (Hel- 
lier et al. 1990; Parmar et al. 2000) and 4U 1700-37 (Rubin 
et al. 1996). 

In the case of non-eclipsing X-ray binaries, such as Cyg 
X-3, the orbital evolution has been measured using its stable 
orbital modulation light curve (Singh et al. 2002). Van der 
Klis et al. (1993) also studied the orbital evolution of the 
LMXB 4U 1820-303 by analyzing the modulation of the 
685 s orbital light curve. 

Amongst the above mentioned LMXBs, the orbital sep- 
aration has been found to be increasing in 4U 1822—37 (Par- 
mar et al. 2000), X2127-I-119 (Homer & Charles 1998) and 
SAX J1808.4-3658 (Jain et al. 2007). The observed orbital 
evolution timescale of ~10® yr in the case of 4U 1822—37 
(Heinz & Nowak 2001) and X2127-M19 (Homer & Charles 
1998) is assumed to be due to short lived mass exchange 
episodes. An orbital evolution timescale of about 70 x 10^ 
yr has been determined in the accretion powered millisecond 
X-ray pulsar SAX J1808.4-3658 (Jain et al. 2007; Burderi 
et al. 2009; Hartman et al. 2009). The timescale of orbital 
evolution was proposed to be due to a strong tidal inter- 
action between the components of the binary. A decreasing 
orbital period has been detected in 4U 1820—303 (van der 
Klis et al. 1993) and Her X-1 (Deeter et al. 1991; Paul et 
al. 2004). However, the observed orbital period evolution is 
faster than that predicted by theoretical models of conser- 
vative mass transfer. In the case of EXO 0748—676, Wolff 
et al. (2009) found four distinct orbital period epochs in the 
last 20 years and attributed it to be due to magnetic cycling 
in the companion star. 

4U 1822-37 is an LMXB with an orbital period of 5.57 
hr (Hcllier et al. 1990; Parmar et al. 2000). It is one of the 



few LMXBs that harbors an accretion powered X-ray pulsar 
(Jonker &: van der Klis 2001). The source is surrounded by 
an accretion disc corona formed by evaporation of matter 
from the inner accretion disc by radiation pressure of the 
neutron star (White & Holt 1982). The light curve exhibits 
a narrow and a broad dip in the intensity. The narrow dip 
in the X-ray light curve is attributed to the partial eclipse of 
the corona by the companion star (White et al. 1981; White 
& Holt 1982; Mason & Cordova 1982; Hellier & Mason 1989, 
Hellier et al. 1992), whereas the broad dip is interpreted as 
resulting from occultation of the corona by a bulge on the 
outer edge of the accretion disc (White & Holt 1982; Hel- 
lier & Mason 1989). Being an eclipsing system, the orbital 
inclination is also known with a small uncertainty (Heinz 
& Nowak 2001; Jonker et al. 2003). An accurate ephemeris 
of 4U 1822—37 is also known from the eclipse timing (Par- 
mar et al. 2000) and the size of the binary orbit is known 
from the pulse timing (Jonker & van der Klis 2001). It is 
thus an ideal system for determination of the masses of the 
stellar components. From the spectroscopic measurements 
of the binary system and by assuming a mass of 1.4 M© for 
the neutron star, Cowley et al. (2003) calculated the mass 
of the companion. A good estimate of the mass of the neu- 
tron star is also known from the K-correction of the radial 
velocity curves (Mufioz— Dariaz et al. 2005). 

We report here on a detailed timing analysis of the 
eclipsing X-ray binary pulsar, 4U 1822—37. We have deter- 
mined 16 new mid-eclipse time measurements of the source 
with data obtained from instruments onboard RXTE, Swift, 
XMM-Newton and Chandra observatories. We have also 
measured the spin period from several newer RXTE obser- 
vations. Combining these new results with previously pub- 
lished measurements, we have derived new updated esti- 
mates of the orbital and spin parameters. Using these mea- 
surements, we have also determined the accretion luminos- 
ity and the magnetic field strength of the pulsar, which is 
crucial for this system. Even though the binary parameters 
have been reliably measured in the past, only rough esti- 
mates of the source intrinsic luminosity and the neutron 
star's magnetic field have been made by various authors, 
with inconsistent results. Assuming a distance of 2 kpc. Ma- 
son & Cordova (1982) had estimated an isotropic luminosity 
of ~ 10^® ergs s^^ . From the broad band spectrum analysis, 
Parmar et al. (2000) measured the 1 — 10 keV X-ray fiux to 
be 5 X 10~^" ergs cm~^ s~^. This implies an X-ray luminos- 
ity of 5.6 X 10^* {d/lkpcf' ergs s~^. Assuming a magnetic 
field strength in the range (1 — 5) x 10^^ G, Jonker & van 
der Klis (2001) derived an intrinsic luminosity of (2—4) x 
10''^ ergs s^^ from the spin-up measurements. 



2 OBSERVATIONS AND ANALYSIS 

In this work, we have analyzed observations made with 
the Proportional Counter Array (PCA) and the All Sky 
Monitor (ASM) on board the Rossi X-ray Timing Explorer 
{RXTE); the X-ray Telescope (XRT) on board Swift; the 
European Photon Imaging Cameras (EPIC) - MOS instru- 
ments on board XMM-Newton; and the Advanced CCD 
Imaging Spectrometer (ACIS) on board the Chandra obser- 
vatory. A log of the X-ray observations used for the present 
work is presented in Table 1. 
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Table 1. Log of observations analyzed in the present work 



Instrument 


Time range 


Observation 


Exposure 




(MJD) 


IDs 


(ks) 


RXTE-ASM 


50088 - 54756 




5096 


RXTE-PCA 


51018 


30060-02-02-00 


18 




51975 


50048-01-01-05 


8 




52094 


50048-01-03-01 


6 




52095 


50048-01-03-02 


8 




52432 


70037-01-01-00 


26 c 




52489 


70037-01-03-00 


21 




52489 


70037-01-03-01 


7 




52491 


70037-01-04-00 


18 




52503 


70037-01-06-02 


4 




52503 


70037-01-06-04 


12 




52519 


70037-01-07-00 


24 




52882 


70037-01-11-00 


11 




52883 


70037-01-11-01 


13 


Chandra-AClS 


51779 - 51780 


671 


40 


XMM-Newton-MOS 


51975 - 51976 


0111230101 


54 


Swift-XKT 


54317 - 54318 


00036691003 


17 ] 



4U 1822-37 was monitored regularly by the RXTE- 
ASM (Levine et al. 1996), which comprises three wide-field 
scanning shadow cameras (SSCs) which are mounted on 
a rotating boom. The SSCs are rotated in a sequence of 
"dwells" with an exposure typically of 90 s, so that most of 
the sky can be covered in one day. The dwell data are also 
averaged for each day to yield a daily-average. The data used 
in the present analysis, covered the time between MJD 50088 
to MJD 54756. The 1.5-12 keV long term ASM light curve 
was corrected for the earth motion using the tool earthSsun 
of the HEAsoft analysis package, ftools ver 6.5.1. The light 
curve, binned with a binsize of 50 days, is shown in Fig- 
ure 1. It shows a gradual decrease in the X-ray intensity 
over the last 13 years, alongwith long term fluctuations over 
timescales of about an year. There could be many reasons 
for the observed decrease in the X-ray flux; decrease in mass 
accretion rate and a change in the structure of the corona 
(Bayless et al. 2010) being two possibilities. 

To check whether the observed decrease in the source 
intensity is an instrumental effect, we also studied the light 
curve of the well known Crab Nebula. As shown in Figure 1, 
the intensity of Crab does not show a decreasing pattern over 
the same time period as in 4U 1822-37. The RXTE-ASM 
measurements of long term intensity variations are quite re- 
liable and for some sources correlated variability has also 
been reported with different instruments (4U 1626-67: Jain 
et al. 2010). 

We have also analyzed the data taken with RXTE- 
PCA, consists of flve xenon/methane proportional counter 
units (PCUs) and is sensitive in the energy range of 2—60 
keV with an effective area of 1300—6500 cm^, depending on 
the number of operating PCUs (Jahoda et al. 1996). Data 
for the eclipse timing analysis were chosen such that it cov- 
ered the entire eclipse phase. The data were taken from the 
Standard- 1 mode of PC A and the background count rate 
estimated using the runpcabackest tool was subtracted from 
the light curves. The photon arrival times of the background 
subtracted light curve were then corrected for the solar sys- 
tem barycenter using the ftool fxbary. To search for the pul- 
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Figure 1. 1.5-12 kcV RXTE-ASM light curve of 4U 1822-37 
and Crab, binned with 50 d. An overall decrease in the intensity 
is clearly seen in 4U 1822-37, whereas the count rate is almost 
constant in case of Crab (Count rate for Crab has been scaled- 
down by a factor of 100). 



sations, we have used the Event mode data with a time 
resolution of 16 ^s. 

4U 1822-37 was observed with XMM-Newton (Jansen 
et al. 2001) for 54 ks on March 07, 2001 (obs ID - 
0111230101). XMM-Newton carries three X-ray mirrors and 
three focal plane instruments, a European Photon Imag- 
ing Camera (EPIC) - pn, MOSl and MOS2, each with a 
field of view of about 30' x 30'. All the cameras (Struder 
et al. 2001; Turner et al. 2001) were operated in full frame 
mode with medium filter. The observation details are sum- 
marized in Table 1. The EPIC observation data files were 
processed using the XMM - Science Analysis System (SAS 
version 8.0.0). For the present analysis, we have used data 
taken with the EPIC-MOS2 in the energy range 0.2-15 
keV. The X-ray events were extracted from a circular re- 
gion of radius 10" centered on the position of the target in 
the EPIC-MOS image. The background X-ray events were 
extracted from a source- free circular region, with a radius of 
20". The background subtracted light curve was barycenter 
corrected using the SAS tool barycen, using the JPL-DE405 
ephemerides. 

4U 1822-37 was observed with Chandra-AClS (Weis- 
skopf et al. 2000) on August 23, 2000 for an exposure time 
covering 2 full binary periods (Cottam et al. 2001). We used 
the Chandra Interactive Analysis of Observations (CIAO) 
software (ver. 4.0; CalDB ver. 3. 4. 2) and standard Chandra 
analysis threads to reduce the data. No background flares 
were found, so all data were used for further analysis. For 
the present work, we used the energy range 0.3—12.0 keV. 
Light curves were extracted from a circular region with a 
radius of 5". Background events were obtained from an an- 
nular region with an inner (outer) radius of 15" (30"). The 
background subtracted light curve was barycenter corrected 
using the CIAO-tool axbary. It should be remarked that pile- 
up is an important phenomena which is inherent to CCD 
detectors, specially with instruments on board Chandra. It 
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is a major concern while measuring fiux/spectra, especially 
when the source is bright. But there will be no effect on the 
determination of the mid-eclipse time. Pile-up can make the 
eclipse shallower but cannot change the mid-eclipse time. 

We have also analyzed data from the Swift observatory 
(Gehrels et al. 2004). The scientific payload consists of a 
wide field instrument, the gamma ray Burst Alert Telescope 
(BAT; Barthelmy et al. 2005) and two co-aligned narrow 
field instruments: the X-ray Telescope (XRT, Burrows et al. 
2005) operating in the 0.2-10 keV energy band and the Ul- 
traviolet/Optical Telescope (UVOT, Homing et al. 2005). 
For the present work, the XRT data was processed with 
XRTDAS software data pipeline package (XRT-PIPELINE 
V. 0.12.0). Calibrated and cleaned level 2 files were produced 
with the xrtpipeline task. We have used an energy range of 
0.2—10 keV and all data were extracted in the Window Tim- 
ing mode, with a total exposure time of 17.3 ks. X-ray events 
from within a rectangular region of width 6 and height 40 
pixels were extracted for timing analysis. Background data 
was extracted from a neighbouring source free region of sim- 
ilar dimensions. We have applied earth2sun correction to the 
background subtracted light curve which was produced with 
a timing resolution of 1 s. 

2.1 Eclipse-timing analysis 

The light curves obtained from Chandra, XMM-Newton, 
Swift and RXTE-PCA observations were folded with the 
known orbital period of 5.5706 hr (Parmar et al. 2000). 
But since the _RXr_E-ASM dataset represents an average 
value over years of data, we folded the 7?Xr_B-ASM light 
curves with the best measurement of the orbital period 
derived from the actual RXTE-ASM data. We obtained 
an orbital period of 20054.27 s and 20054.24 s from the 
RXTE-ASM data spanning MJD 50088-52432 and MJD 
52432-54756, respectively. The RXTE-ASM light curve was 
therefore, folded in two segments (MJD 50088-52432 and 
MJD 52432-54756), while a single folded profile was ob- 
tained from each of the other observations. 

The folded light curves of the data obtained from 
RXTE-ASM, Chandra, XMM-Newton and Swift missions 
are shown in Figure 2. The RXTE-PCA light curves are 
shown in Figure 3. It should be noticed that the eclipse mor- 
phology is changing. The eclipses seem to vary in depth and 
shape. The variable eclipse depth shows that the projected 
geometry of the accretion disk and corona is changing. It is 
also possible that a part of the change in the eclipse mor- 
phology in Figure 2, is due to the slightly different energy 
bands used and the differences in instrument efficiency over 
the energy bands. 

The light curves show clear signs of orbital modulation 
(i.e. partial eclipse and a sinusoidal modulation) with an 
orbital period of 5.57 hr. A model consisting of a Gaussian 
and a constant was fit to the eclipse interval (0.45—0.55 
orbital phase) in each folded light curve (as in Parmar et al. 
2000). The reduced x^ of the fits were in the range 0.5—4 
for 13 d.o.f. Figure 4 shows the eclipse interval (0.45—0.55 
orbital phase) for one of the folded RXTE-PCA light curves 
(Obsid 70037-01-03-00). Gaussian and a constant model was 
fit to this interval. Solid line in the top panel of Figure 4 
shows the best fit model and the bottom panel shows the 
residuals of the fit. 
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Figure 2. Folded light curves of 4U 1822-37 from the data ob- 
tained from RXTE-ASM, Chandra, XMM-Newton and Swift- 
XRT. The Chandra, XMM-Newton and Swift-XKT light curves 
were folded with a period of 20054.3 s, while the RXTE-ASM 
light curves were folded with 20054.27 s (MJD 50088-52432) and 
20054.24 s (MJD 52432-54756). The epoch of the mid eclipse time 
for each observation is given in each panel (in MJD). The Y-axis 
shows the intensity normalized by the average intensity and on 
the X-axis, two orbital cycles are shown for clarity. 



The arrival time of the eclipse which occurred closest to 
the mid time of the observation was taken for further anal- 
ysis. The new mid eclipse time measurements along with la 
uncertainties are given in Table 2. The orbit number (cy- 
cle) is with respect to the first reported mid-eclipse time 
(Hellier & Smale 1994). The newly determined arrival times 
were combined with the earlier known values and fitted with 
a quadratic model. We obtained a x^ of 432 for 35 d.o.f. 
However, it should be remarked that the uncertainties in 
the newly determined mid eclipse time measurements are 
too small to give a reliable estimate. Therefore, we rescaled 
the errors in the individual measurements (by multiplying 
the individual errors by the square root of the above men- 
tioned reduced x^); in order to compare the results with 
the earlier known estimates of the orbital parameters (Hel- 
lier & Smale 1994; Parmar et al. 2000). The 16 new arrival 
times were combined with the earlier observations and were 
again fitted with a quadratic model to obtain an updated 
ephemeris. The best fit model gave an orbital period (Porb) 
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Figure 3. 2-60 keV RXTE-PCA light curves of 4U 1822-37, 
folded with a period of 20054.3 s. The epoch (in MJD) of each 
observation is given in the corresponding panel. The Y-axis shows 
the intensity normalized by the average intensity and on the X- 
axis, two orbital cycles are shown for clarity. 



of 0.23210887(15) d and a period derivative (Port) of 1.3(3) 
X 10""' d d"^ (at To = MJD 45614) with a x^ of 36.2 for 35 
d.o.f. This imply an orbital evolution timescale of 4.9(1.1) x 
10^ yr. The resulting best fit orbital parameters are given in 
Table 3. We subtracted the best fit linear component from 
the ephemeris history and the residuals are plotted in Figure 
5. The derived values of the rate of change in orbital period 
and the timescale of orbital evolution are consistent, within 
measurement errors, with those obtained from timing of the 
eclipses in the optical and UV data (Bayless et al. 2010). 



2.2 Pulse-timing analysis 

We have performed a pulsation analysis to determine the 
spin period of the neutron star and the pulse period evo- 
lution. The light curves were corrected for the orbital mo- 



Table 2. New X-ray mid-ecUpse times of 4U 1822-37. 



Cycle MJD 



Uncertainty (d) 



Satellite 



23281 


51018.54400 


0.00006 


RXTE-PCA 


24300 


51255.0638 


0.0018 


RXTE-ASM-\ 


26562 


51780.0954 


0.0004 


Chandra 


27406 


51975.9968 


0.0001 


RXTE-PCA 


27407 


51976.22941 


0.00004 


XMM-Newton 


27918 


52094.83802 


0.00007 


RXTE-PCA 


27922 


52095.76593 


0.00006 


RXTE-PCA 


29374 


52432.78961 


0.00005 


RXTE-PCA 


29618 


52489.42192 


0.00005 


RXTE-PCA 


29628 


52491.74614 


0.00007 


RXTE-PCA 


29678 


52503.34773 


0.00005 


RXTE-PCA 


29747 


52519.36724 


0.00006 


RXTE-PCA 


31311 


52882.38623 


0.00006 


RXTE-PCA 


31315 


52883.3151 


0.0001 


RXTE-PCA 


34378 


53594.2680 


0.0020 


RXTE-ASM-\ 


37498 


54318.4506 


0.0004 


Swift-XRT 



fTlie mid eclipse times represent the time span MJD 
50088-52432 and 52432-54756. 
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Figure 4. The expanded view of the eclipse interval (0.45—0.55 
orbital phase) of one of the folded RXTE-PCA light curve of 4U 
1822-37 (Obsid 70037-01-03-00; Figure3(f)). A model consisting 
of a Gaussian plus a constant was fit to this interval. The solid 
line in the top panel shows the best fit model (with a reduced x^ 
of 1.2 for 13 d.o.f.) and the bottom panel shows the residuals of 
the fit. 
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Table 3. Orbital parameters of 4U 1822-37. The numbers in the bracket indicate the la uncertainty. 



Parameter Hellier & Smalc Parmar ei al. 

(1994) (2000) 



Best fit value from 
the present analysis 



To (MJD) 45614.80962(15) 45614.80964(15) 45614.80951(12) 

Forb (d) 0.232108788(64) 0.232108785(50) 0.23210887(15) 

c (10-" days)t 2.37 ± 0.56 2.06(23) 1.53(40) 

^orb/Porb (106 yr) 3.1 ± 0.7 3.6 4.9(1.1) 



tQuadratic Ephemeris; JD = Tq + NxPo,.;, + N^ X c 

Table 4. Mid-eclipse times and the corresponding pulse period. 
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Figure 5. The mid eclipse time residuals of 4U 1822-37 have been 
plotted as a function of the orbital cycle, relative to the best fit 
linear ephemeris (orbital period of 0.23210854 d). The solid curve 
is the best fit of a quadratic function having a x^ of 36.2 for 35 
d.o.f. The eclipse times are tabulated in Table 2. The curvature 
of the locus of the residuals is a measure of the orbital period 
derivative of the binary system. The square boxes are the values 
known from previous measurements. The newly determined mid 
eclipse times are shown with "o" . The two horizontal bars indicate 
the time span of the RXTE-AS,M data. 



Mid-eclipse times (MJD) Spin period (s) 



51975.9968(1) 
52094.83802(7) 
52095.76593(6) 
52432.78961(5) 
52489.42192(5) 
52491.74614(7) 
52503.34773(5) 
52519.36724(6) 
52882.38623(6) 

52883.3151(1) 



0.59290132(11) 
0.59286109(8) 
0.59286421(12) 
0.5927922(13) 
0.5927790(06) 
0.5927795(11) 
0.5927737(10) 
0.5927721(08) 
0.5926793(15) 
0.5926852(21) 



completeness. A considerable variation was seen in the pulse 
profile. The pulse profile exhibited a non-sinusoidal profile 
at To = MJD 51975; a relatively broad maximum at Tq = 
MJD 52432; a sharper profile at around To = MJD 52491 
and a sinusoidal variation around To = MJD 52519. How- 
ever, it is difficult to quantify the observed variation in the 
pulse shape at this point. 

We also tried to obtain an independent measurement of 
orbital evolution using the technique of pulse folding and x^ 
maximization (Nailc & Paul 2004; Jain et al. 2007). But it 
should be remarked that in this source, the light travel time 
across the orbit {a^ sin i) is only a factor of two larger than 
the spin period. Moreover, in case of pulse timing technique, 
the pulse profile is assumed to be invariant. But even a small 
orbital phase dependence of the pulse shape, for example 
caused by varying absorption can lead to systematic errors 
in measurement of the orbital parameters. Therefore, though 
this analysis has proved to be successful for other LMXBs, 
we could not find a very accurate measurement of the orbital 
parameters by this technique. 



tion using the long term orbital solution obtained from the 
eclipse timing technique described above. Figure 6 shows 
the spin period history of the neutron star. The pulse pe- 
riod was found to be continuously decreasing with time at an 
average rate (Pspin) of -2.481(4) x 10~^^ s s~^, indicating 
a spin-up timescale of 7578(13) yr. From the observations 
made in 1996 and 1998, Jonker & van der Klis (2001) had 
calculated an average spin-up {Pspin/^spin) of -1.52 x lO^** 
yr^^. The pulse profiles obtained after correcting for the or- 
bital motion, were normalized by the average intensity and 
are shown in Figure 7. The value of the spin period, deter- 
mined by Jonker & van der Klis (2001) are also marked for 



3 DISCUSSIONS 

We have performed a detailed timing analysis of the low 
mass X-ray binary pulsar 4U 1822—37 using the archival 
X-ray data from several X-ray observatories. Observations 
used in the present work covered a time span of 13 years 
and more than 14,000 binary orbits. Using the 16 new ac- 
curately measured mid-eclipse times, we have obtained an 
orbital period of 0.23210887(15) d with a significant orbital 
period derivative of 1.3(3) x 10"^" d d"^ (at To = MJD 
45614). It indicates an orbital evolution timescale Vorb/Porb 
= 4.9(1.1) Myr. The orbital and spin parameters were also 
measured by correcting the light curves for the binary mo- 
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Figure 6. The variation of spin period (Pspin) of the neutron 
star with time. The first two points indicate the spin period de- 
termined by Jonker & van der Khs (2001). The straight hne rep- 
resents the best fitted Unear curve, with a x^ of 13.7 for 12 d.o.f. 



tion of the pulsar and then optimizing the pulse detection. 
However, the results from the pulsation analysis did not im- 
prove the orbital evolution measurements. 



3.1 Orbital Evolution 

X-ray binaries can evolve by various mechanisms such as 
mass transfer within the system due to Roche lobe overflow, 
tidal interaction between the components of the binary sys- 
tem, gravitational wave radiation, magnetic braking, and 
X-ray irradiated wind outflow (As mentioned before in Sec- 
tion 1). Orbital evolution has been measured in some other 
low magnetic field LMXBs, such as 4U 1820—30 (van der 
Klis et al. 1993), EXO 0748-676 (Wolff et al. 2008) and 
SAX J1808.4-3658 (Jain et al. 2007), and several models 
have been proposed to explain the orbital period evolution 
in them. In the case of 4U 1822-37, the reason behind a high 
rate of orbital evolution is not known. The measured rate of 
change of orbital period in 4U 1822—37 is much greater than 
that expected due to gravitational wave radiation (Verbunt 
1993). We also note that the timescales of orbital evolution 
due to tidal interaction between the components of the bi- 
nary system ranges from a few Myr in HMXBs to about 
10^° years in LMXBs (Applegate & Shaham 1994). Conser- 
vative mass transfer, mass loss from the binary system due 
to an X-ray irradiated wind outflow and magnetic cycling in 
one of the binary components are the other possibilities. We 
first examine the possibility and consequences of the case 
of a conservative mass transfer in this system, followed by 
the case of an X-ray irradiated wind outflow and magnetic 
cycling in the companion star. 
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Figure 7. Normahzed X-ray pulse profiles of 4U 1822-37, ob- 
tained from observations with RXTE-PCA. We have divided the 
pulse period into 32 bins. The epoch (in MJD) of each observation 
is mentioned in the corresponding panel. 



3.1.1 Conservative mass transfer 

In the case of conservative mass transfer, the mass transfer 
rate from the companion star equals the accretion rate onto 
the neutron star. We have estimated the mass accretion rate 
using the best known estimate of the mass of the neutron 
star (M„s) and the companion star (Mc) (Mufioz-Dariaz et 
al. 2005) in this system: 



1.61Mq < M„s ^ 2.32Mq 
O.44M0 < Mc ^ 0.56Mq 

and a mass ratio q ( 

0.24 s; g s; 0.27 






in the range: 



(1) 
(2) 

(3) 



In the case of compact binaries, the orbital angular momen- 
tum (J) is given by (King 1988): 



J = McM,^, 



where, 



^Gay 
km) 



(4) 



M = Total mass of the binary system 
G = Gravitational constant 
a = Binary orbital separation. 



Mc + M„ 
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For conservative mass transfer, all the mass lost from 
the companion is accreted by the compact object (M„s + 
Mc = 0) and the orbital angular momentum is conserved 
(J = 0). In such a case, the mass accretion rate {Mns) is 
related to the orbital period [Porb) as: 



Mr,s = 



3 "Porb \ Mns) 



(5) 



Since the observed -^^ is 2.10 x 10 ^ yr ^, therefore, 

^orh 

we get a mass accretion rate of (4.2—5.2) x 10^* M0 yr^^ 
and a mass transfer timescale of ~ 10 yr. Assuming an 
accretion efficiency {if) of 10%, this implies an accretion lu- 
minosity hacc {= rime's/2) in the range (1.19-1.48) x 10** 
ergs s~ . 

For a neutron star with a magnetic moment /i (~ Br"^), 
the spin frequency derivative v is related with the mass ac- 
cretion rate as (Frank et al. 2002): 



z> ~ 0.47 X 10 
where. 



-12 3/7 



m,' mi6 



6/7 2/7 J--1 



2/7r-lTT -1 

M30 hr, Hz s 



(6) 



mi — mass (m) of the neutron star in units of 1 Mq 

niiQ— accretion rate in units of 10^^ g s^^ 

/iao ~ magnetic moment in units of lO"'" G cm* 

I45 — the moment of inertia of the accreting star in units of 

10*^ g cm^. I45 ~1 for a neutron star. 

Using data obtained from RXTE-PCA, we have de- 
tected 0.59290132(11) s (at To = MJD 51975) X-ray pul- 
sations with an average spin-up rate of -2.481 x 10^^^ s 
s^^. These values imply a magnetic field strength of (1—3) 
X 10* G. 

The following caveats apply to the case of conservative 
mass transfer: 

• It requires a large mass transfer which corresponds to 
luminosity near the Eddington rate. Considering the un- 
certainties in some of the parameters used above, a near- 
Eddington accretion luminosity is not inconceivable, espe- 
cially because this source has a corona surrounding it (White 
& Holt 1982; Heinz & Nowak 2001; Bayless et al. 2010). A 
comparison of the L^/ljopt ratio of this system with other 
LMXBs also suggest that the true X-ray luminosity of the 
central X-ray source is probably significantly higher. See 
Bayless et al. (2010) for more discussion on this. For exam- 
ple, in SS433 the mass transfer rate is believed to be much 
higher than the Eddington rate, that results in accumula- 
tion of material around the compact object which blocks 
its X-ray emission (Begelman, King, & Pringle 2006; Clark, 
Barnes, & Charles, 2007). 

• It requires the neutron star to have a low magnetic field 
strength of (1—3) x 10* G while from the high energy cut- 
off in the X-ray spectrum (Parmar et al. 2000), Jonker &: 
van der Klis (2001) determined a magnetic field strength of 
~(1— 5)xlO^'^G. It should be noted that the coronal X-rays 
dominate the X-ray spectrum of this source and the pulsed 
X-rays contribute to only a few percent of the total X-ray 
emission. Thus the X-ray spectral shape is unlikely to be a 
reliable indicator of the magnetic field strength of the com- 
pact star. 

• In accreting neutron stars, we expect pulsations only if 
the compact object possesses a magnetic field strong enough 
to disrupt the inner regions of the accretion disc and chan- 



nel the accretion flow onto the polar caps. In other words, 
accretion onto the neutron star is controlled by the mag- 
netic field if the magnetospheric radius is larger than the 
stellar radius. In case of 4U 1822-37, the observed values of 
the mass transfer rate and the magnetic moment, indicate a 
magnetospheric radius of 2 x 10® cm, which is of the same 
order as the neutron star radius. Therefore, it is not certain 
whether the conservative mass transfer is indeed responsi- 
ble for the observed changes in the orbital period. However, 
several low mass X-ray binaries, such as SAX J1808. 4-3658, 
do show pulsations at similar luminosity level even though 
they have a low surface magnetic field of (1—5) x 10* G (Di 
Salvo & Burderi 2003). 

• For the accretion luminosity, spin period and mag- 
netic field strength mentioned above for a conservative mass 
transfer case, the neutron star is far from spin equilibrium. 
Detection of such a system is a novelty and unlikely, un- 
less the current X-ray state is a long lived transient phase. 
Here one may note that in the recent years, many LMXBs 
have been found which spend a small fraction of the time in 
transient high states, for example, the millisecond accreting 
pulsars. 



3.1.2 X-ray irradiated wind outflow 

Mass loss from the binary system can occur in the form of 
X-ray irradiated wind outflow. In case of LMXB 4U 1822-37, 
there is no signature of wind outflow in the form of absorp- 
tion lines in high resolution X-ray spectrum (Gottam et al. 
2001). Recently, Bayless et al. (2010) reported a broad C IV 
emission line in the UV spectrum of 4U 1822-37, indicating 
a strong disk wind outflow. They have derived a wind out- 
flow velocity of 4000 km s^^, based on a measured width of 
45 A° of the C IV emission line. However, these observations 
arc insufficient to estimate the total mass outflow rate. 



3.1.3 Magnetic cycling in the Companion star 

Secular changes such as magnetic cycles in the secondary 
star (Hellier et al. 1990) is also a possible mechanism re- 
sponsible for the observed orbital evolution in 4U 1822-37. 
However, the spectral type and the evolutionary history of 
the companion star is unknown (Muiioz- Dariaz et al. 2005). 



4 CONCLUSIONS 

We have presented new and more accurate measurement of 
orbital evolution of the LMXB 4U 1822-37 and a longer time 
base for measurement of its spin evolution. Considering the 
possibility of a large intrinsic L^, and some evidence of mass 
outflow, we conclude that the orbital evolution in this sys- 
tem is complex including the effects of a large mass transfer 
rate and X-ray irradiated wind outflow. In this scenario, the 
magnetic field strength of the neutron star is probably in 
between the same in typical low mass X-ray binaries and 
high mass X-ray binaries. 
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